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Based on the self-consistent scheme beyond mean-field approximation in 
the large N c expansion, including current quark mass explicitly, a general 
scheme of SU(2) NJL model is developed. To ensure the quark self-energy ex- 
panded in the proper order of N c , an approximate internal meson propagator 
is deduced, which is in order of 0(1/N C ). In our scheme, adopting the method 
of external momentum expansion, all the Feynman diagrams are calculated 
in a unified way by only expanding the quark propagator. Our numerical 
results show, that different from the mean field approximation in which the 
explicitly chiral symmetry breaking is invisible, the effect of finite pion mass 
can be seen clearly when beyond mean-field approximation. 
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I. INTRODUCTION 



As we all know, chiral symmetry breaking was originally explained quite well by Nambu- 
Jona-Lasinio (NJL) model as early as 1961 0. Like in superconductivity, the strong attrac- 
tive force between quark and antiquark in the J p = + channel can create non-perturbative 
ground state with qq condensation. Due to the pair condensation, the original symmetry be- 
tween massless left and right-handed quarks is broken down to Uv(Nf), and then the quarks 
obtain constituent mass. The remaining residual interactions between the constituent quarks 
bind them into collective excitations, i.e., hadrons in the chiral symmetry breaking vacuum. 
Especially in the pseudoscalar channel the residual strong interaction creates massless pions 
as Goldstone bosons in the chiral limit. When a small current quark mass mo is introduced 
in the theory, chiral symmetry is explicitly broken, and pion obtains its small physical mass 
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Although the NJL model has two serious drawbacks, i.e., lacks of confinement and renor- 
malizability, it is still regarded as an applicable model at low momentum, especially for 
dealing with processes of pion, such as pion-pion scattering near threshold. 

Traditionally, the scheme of the NJL model is represented by two Schwinger- Dyson (SD) 
equations, one is for the constituent quark propagator, and the other is for the composite 
meson propagator. At the lowest level, the applications of the NJL model are based upon 
mean-field approximation @ - @, i.e., Hartree approximation to the gap equation for quark 
mass and the random-phase approximation (RPA) to the Bethe-Salpeter equation for meson 
mass. It is clear, that at this level the solution of the gap equation determines the meson 
propagators, but the solution of meson SD equation has no feedback to the quark propa- 
gator. Since the constituent quark propagator is the fundamental element, from which all 
the quantities, including quark mass, meson masses and quark-antiquark condensate, are 
calculated, it is necessary to consider the back contribution of meson modes to the quark 
propagator. 



Among efforts [Bf - WM to go beyond the mean-field approximation, Refs. [O and (13 
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are in a chirally symmetric self-consistent approximation, namely the chiral properties such 
as the Goldstone's theorem, the Goldberger-Treiman relation and the conservation of the 
quark axial current are exactly preserved in the chiral limit of the NJL model. By using 
effective action method in a semi-bosonized way, and expanding the action to one quark- 
loop and one- meson- loop in ||12||, or directly evaluating the Feynman diagrams under the 



constraint to keep the chiral relations at quark level in |13 . 

In this paper, we extend the method of |13[ to a general scheme with explicit chiral 
symmetry breaking in the SU(2) NJL model. Different from the case in the chiral limit, 
we must be careful to deal with the form of internal meson propagators. In a way different 
from fl3 |, we regard the constituent quark as the fundamental element and only expand 



quark's propagator in the power of small external momentum in the calculation of Feynman 
diagrams. 

In the process to go beyond the mean-field approximation, we have to deal with the 
divergent integrals of quark loops and meson loops. We adopt Pauli-Villars regulation PJT5 



to treat divergent integrals resulted from quark loops, and choose a covariant cutoff A& for 
the meson momentum. There are four parameters in our treatment, namely the current 
quark mass m , quark coupling constant G, fermionic cut-off Aj and bosonic cut-off A b , 
to be fixed. In the mean-field approximation, the three parameters mo, G, A/ are usually 
fixed by comparing with the pion mass = 140 MeV, pion decay constant f n = 92.4 
MeV and the quark condensate 1/2 < qq > 1 ' 3 = —250 MeV. In the near future, the DIRAC 
experiment will measure the ir — ir scattering lengths in good precision, which will shed some 
light on the quark condensate |L6|]. To see clearly the quark condensate dependence of the 
four parameters, we give only the quark condensate a reasonable constraint: -300 MeV ~ 
-200 MeV. 

The outline of this paper is as follows: In section 2, we briefly review the general scheme 
represented by two Schwinger- Dyson equations in the SU(2) NJL model. In Section 3, 
we introduce the method of external momentum expansion, and prove a general relation 
between the pion polarization function and the axial- vector matrix element. We also deduce 



the internal meson propagator to 0(1/N C ) order in the N c expansion. Our numerical results 
with mesonic contributions and the effect of explicit chiral symmetry breaking will be shown 
in section 4. The conclusions are given at the end. 

II. THE FUNDAMENTAL SCHEME IN SU(2) NJL MODEL 

A. Two Schwinger-Dyson equations in SU(2) NJL model 

In this section, we briefly review the traditional scheme of SU(2) NJL model with large 
N c expansion. The two-flavor NJL model is defined through the Lagrangian density, 

C = ${i-fd» - m )^ + G[{^f + (^75T^) 2 ], (1) 

here G is the effective coupling constant of dimension GeV -2 , and m is the current quark 
mass, assuming isospin degeneracy of the u and d quarks, and i/j,iJj are quark fields with 
flavor, colour and spinor indices suppressed. 

The traditional non-perturbative method of NJL model is inspired from many-body 
theory. The complete description is represented by two Schwinger-Dyson (SD) integral 
equations, i.e., the constituent quark propagator, see Fig. la, and the composite meson 
propagator, see Fig. lb. 
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FIG. 1. The two Schwinger-Dyson equations for the quark propagator a and the meon propa- 
gator b. The light dashed lines in a and b represent the four-fermion vertex 2iG, and the one-vertex 
gray bubble in a and the two-vertex gray bubble in b indicate the quark self-energy and the meson 
polarization function respectively. 
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The two SD equations must couple to each other self-consistently and keep chiral symmetry 
in the chiral limit. The one- vertex grey bubble kernel a in Fig. la represents the quark self- 
energy, and the two-vertex grey bubble kernel b in Fig. 16 indicates the meson polarization 
function. In the case of SU(2), meson modes refer to only pseudoscalar mesons tt and scalar 
mesons a. 

It is difficult to give the full expressions of kernel a and kernel b in Fig. 1. Usually 
an approximation scheme called large N c expansion is adopted in NJL model, i.e., the both 
kernels expanded in 1/N C , where N c is the number of color. In this scheme, GN C is a 
constant when N c — > oo, so coupling constant scales like G ~ N~ x . One can find the 
detailed description of large N c expansion in |T7|,|9|. V.Dmitrasinovic et.al. proved in their 



paper [|1^] that the kernel a and kernel b shown in Fig. 2 are self-consistent leading and 
subleading order in N c expansion and can keep all the chiral relations in the chiral limit. 

The leading 0(1) order of quark self-energy kernel a in Fig. 2 is named m# after 
Hartree approximation, and the subleading 0(1/N C ) order is called 5m. Including current 
quark mass m , the gap equation for quarks can be expressed as 

m = mo + m>H + 5m, (2) 

where = !6imGN c F, and F is given in Appendix in order to compactly discuss physics 
in the text. 

The leading 0(N C ) and subleading 0(1) order of meson polarization function kernel b 
are expressed as H^ P (k) and oTIm(^) respectively, here M represents tt or a. It is clear 
to see that the back interaction which conserves all the chiral properties is reflected in the 
contribution from the meson propagator to the quark mass. The total meson propagator 
can be expressed in terms of the total polarization function 

where Tl M (k) = n^ PA) (A;) +5U M (k), and the RPA part is Ui RPA \k) = 8iN c F -4iN c k 2 I(k) 
for tt and n^ Pj4) (A;) = 8iN c F - AiN c (k 2 - Am 2 )I(k) for a, and the function I(k) is defined 
in Appendix. 
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FIG. 2. kernel a and kernel b in the quark and meson propagators, mu and Sm are the 
leading and subleading contributions to the quark mass. 

and 5U^ d) are the leading and 

subleading order contributions to meson polarization function. The heavy solid lines indicate the 
constituent quark propagator, and the heavy dashed lines represent n or a propagator —iD^ PA \q) 
in RPA approximation. 



In the NJL model, mesons are bound states of constituent quark and anti-quark, the meson 
mass m M satisfies the meson propagator's pole condition 

I - 2GU M (k 2 = m 2 M ) = 0, (4) 

and the coupling constant guqq is determined by the residue at the pole 

9M 2 qq = {dIi M {k)/dk 2 )-\, =ml . (5) 

Performing direct calculation of the Feynman diagrams in Fig. 2 leads to the explicit 
expressions for the quark mass 

-i5m=12GN c J ^i\-^r A \q)] J ^-Jr[ l5 S(p + q h 5 S(p)lS( P )] 

~ AGNc S J^ [ - iD ° RPA){q)] I ^niS( P + q)lS(p)lS(p)}, (6) 
and for pion polarization fluctuations 



Sn®(k)=2iN e [f J^[-^i RPA) (q)] J -^Tr[^S(p + qh 5 S(p)^S(p-k)^S(p-k + q)} 

+ / 0jl [ - iD ° RPA){q)] I ^n^S(p + q )lS(ph 5 S(p-k)lS(p-k + q )}}, (7) 
5Yl^(k)=-mN c [J 0- 4 [-iDi RPA \q)] J ^Tr[^S(p)^S(p + q h 5 S(p)^S(p-k)} 

~l I J^ [ - iDi * PA){q)] I J^nbtStolSip + q)lS(p) l5 S(p - k)}}, (8) 

^ d) W=i/ ^[-^(Ml^ (9) 

Here, the expression T an7T (k, q) is the (cr7T7r)-vertex 

- ir™(fc, q) = -4N c t J 0-Jr[ l5 S(p + k)lS(p + qh,S(p)], (10) 

S(p) = l/(j6 — m) the quark propagator and — \D^ A \q) the internal meson propagator. 
In all the formulae of this paper, k means the n external momentum, p the internal quark 
momentum and q the internal meson momentum; 



B. Pion Decay Constant 

Another important quantity is the pion decay constant f n which is calculated from the 
vacuum to one-pion axial- vector matrix element, replacing one vertex 275 r of kernel b in 
Fig. 2 by ig nqq (k)j^ fl T/2. Similar to the expression of U n (k), we express f n (k) as 

U(k) = fl RPA \k) + 5U(k) 

= fi RPA \k) + Sfi b \k) + 5jf\k) + 5fi d \k). (11) 

Since the Feynman diagrams for Tl n (k) and f n have the same structure, we expect to 
find a general relation between them. Starting from the axial-vector matrix element 
< 0| J^ 5 |7r J (/c) >= ifnityk^i, and using the single-quark axial Ward identity involving 
the quark propagator S(p), 

7 5 ^f = 2m 75 + l5 S-\p + k) + 5- 1 (p) 75 , (12) 
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one can deduce the general relations functions 



The detailed proof can be found in Ref. ||13|| . Using pion's pole condition Eq. (f|) and gap 
equation Eq. (TJ), we get finally the relation 

m lU = mo ^ 

Qirqq 2G 

where we have defined the on-shell quantities f n = f^ik 2 = m 2 ), g nqq = g^^k 2 = m 2 ). 
While in the chiral limit, /„- satisfies the Goldberger-Treiman relation f^g^qq = m. 



III. EXTERNAL MOMENTUM EXPANSION 

A. Method of External Momentum Expansion 

We are familiar with one-loop Feynman diagrams, but it is difficult to analytically 
calculate the two-loop Feynman diagrams such as those in kernel b. Since pion's external 
momentum k 2 = m 2 is small, we will calculate the subleading corrections in kernel b to k 2 
term only, but keep the complete contributions from the leading order. 

In the NJL model, the constituent quarks are the fundamental elements, and the mesons 
are their bound. In fact, all the quantities in the NJL model are calculated from the quark 
propagator S(p) = — m). Therefor, if we want to build up a self-consistent expansion 
in the caculation of the the 2-loop diagrams, we should only expand the quark propagator 
in the external momentum k 

S (P T fc ) = jfzpfS-m = S{p) ± S W S b) + S(PWS(PWS(P) + ■■■■ (15) 

With this expansion, all the two-loop diagrams can be expanded naturally in k in a unified 
way. 

The Lorentz covariance gives the expansion form of H n (k) as 
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su^k) = su^o) + k 2 sn^(o) + k 4 su^(o) + ■■■. (ie) 

In the following, we consider 5U n (k) to the k 2 term only. In this approximation, the pole 
condition for the pion propagator is simplified as 

1 - 2G{Yi[ RPA \k 2 = ml) + 511,(0) + m 2 jrt 2 \0)) = 0. (17) 

With the expression of H ( ^ iPA \k) and the gap equation Eq. (@) we obtain the gap equation 
to determine the pion mass, 

ml = — m . (18) 

Gm(-8^/K) + 25ITi 2) (0)) 

Similarly, the expansion of the residue of the pole to the order k 2 gives the pion-quark 

coupling constant 

9~l = (dUi RPA \k)/dk 2 )-\ 2=ml + OT^O), (19) 

where (dYl^ RPA \k) / dk 2 )' 1 = -2JV c i(/(0) + I(k) - k 2 K(k)), and the quark-loop integral 
K(k) is defined in Appendix. 

B. Calculation of 8U {2) (0) 

It is tedious to calculate 511^(0) directly. Here we will use another equivalent way 
by calculating the first term of the axial-vector matrix element in the external momentum 
expansion. 

We have shown the general relations between pion polarization function and axial- vector 
matrix element in Eq. fll3|) . If k 2 = 0, i.e., in the chiral limit, Eq. fll3|) becomes 

^ = ^(0), 6 ^ = m6U n (0). (20) 
Thus the Eq. ( |13|) in general case can be rewritten as 



9irqq{ K ) 9nqq{ K ) 
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With the expansion form of U n (k), the second equation of ([H]) becomes 

kHf ^ = m(k 2 SU (2 \0) + k 4 6U^(0) + ■•■)• (22) 

By replacing in the Feynman diagrams in kernel b one vertex ij^T by ig nqq (k)'~f5'j f j,T/2, the 
obtained subleading axial- vector matrix element has the following form in external momen- 
tum expansion: 

MI^ = kfl m(M(0) + k 2 M^(0) + •••)• (23) 



Comparing the two equations (P2|) and (|23|) , we derive the pion polarization fluctuation in 



terms of the first coefficient of the expansion of the axial- vector matrix element. 

SIl^ (0) = M (0) = M (fe) (0) + M (c) (0) + M (d) (0) (24) 

with 

M«(0) = i N c { J ^[-iDi RPA \q)](-3q 2 L(q)) 

+ / ^[-^ RPA \q)](AK(q) + 3(4m 2 - g 2 )7(g))}, (25) 
M (c) (0)= z iV c {6 J -0- A {-iDi RPA \q)](K(q) + 3K(0) - 3q 2 M(q)) 

+2 / ^i-iDl RPA \q)](5K(q) + 3K(0) - 3(q 2 - 4m 2 )M(g))}, (26) 

M«(0)= -32ziV c / ^iV c [-iD^)(g)][-iD(«^)(g)]{-(/(g) + 2m 2 ^(0))(7(g) -7(0)) 
+ (7(g) + 7(0) - {q 2 + 2m 2 )K(g))7(g) - g 2 7(g)(7(g) + 2m 2 K (0))[-iD^ (q)}} , (27) 

where 7(g), M(q) are defined in Appendix. In calculating M^(0), we have expanded the 
internal cr propagator as 

\-iDl RPA \q - h)\ = [-iDi RPA \q)]{l + q^[-W^ (q)][-Wi RPA \q)]\ (28) 

with [-iDW(q)] = 8N c (I(q) + ((Am 2 - q 2 ) /2q 2 )(I(q) - 7(0) + q 2 K(q)) [-iD^ (q)}). 
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C. Internal RPA Meson Propagator 



Till now, we have assumed that kernel a and kernel b shown in Fig. 2 are expanded 
properly to leading and subleading order in N c . In the chiral limit beyond mean-field ap- 
proximation, the leading meson polarization function U^ PA ^ has the order of 0(N C ), which 
gives the RPA meson propagator D < y PA \q) in the order of 0(1/N C ), thus 5m is the sub- 
leading 0(1/N C ) order of quark self-energy. However, it is not the case when current quark 
mass is introduced beyond mean-field approximation. 

Now we first analyze the one-quark- loop pion propagator D^ PA \q) beyond mean- field 
approximation and including m explicitly. The definition of RPA pion propagator is 

~ in ^- 2 cn - (29) 

Substituting the expression of U^ PA ^(q), and using the gap equation Eq. (||), the denomi- 
nator is 

1 - 2GUi RPA \q) = m » + 6m + i 8GN c q 2 I(q). (30) 

m 

With the pion pole condition Eq. (|J), we can deduce 

m ° + 6m = _ z 8GN c m 2 J(m n ) + 2G5U 7T (m 7T ). (31) 
So the RPA pion propagator fl2D|) can be expressed as 

- xDf PA \q) = 4jVc( _ m 2 /(m j + g2/(g)) _ iSU ^ mn ) ■ ^ 
Here 5n ?r (m 7r ) is the subleading diagrams' contribution, which is in the order of 0(1). This 
term induces the RPA propagator expanded to 0(1/A") order, where n could go to oo. If 
we substitute this complete form of internal pion propagator into kernel a, it is clear that 
the 0(1/A C ) part of RPA propagator corresponds to the 0(1/N C ) part of 5m, and other 
suppressed parts of RPA propagator correspond to the suppressed parts of 5m. 

In order to keep the proper subleading order 0(1/N C ) of 5m, the internal meson prop- 
agator must be in its leading order 0(1/N C ). We neglect the contributions resulted from 
5U n (m n ), and get the leading order of pion propagator in 0(l/iV c ), 
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iDl RPA \q) 4]Vc( _ m 2 /(m ^ + q 2 I{q) y ( 33 ) 
Similarly, the 0(1/ N c ) sigma propagator is 

- LDf PA ) (g) = 4Arc( _ m 2j( m7r ) + ( ?2 _ 4m 2 )/(g)) ' 

It can be seen that G is not included explicitly in the two RPA meson propagators, and this 
simplifies our numerical calculations. 

IV. NUMERICAL RESULTS AND DISCUSSION 

Now we turn to the numerical evaluation. As evaluated in the last two chapters, we 
have three equations, the gap equation, Eq. (g), pion pole condition, Eq. (f|), and pion 
decay constant, Eq. (0). I the chiral limit, pion pole condition Eq. (f|) is trivial and the 



pion decay constant satisfies the Goldberger-Treiman relation. 

In order to deal with the divergence resulted from quark-loop integral and meson-loop 
integral, like in [[R|], we introduce two cut-off, the quark momentum cut-off Af in Pauli- 
Villars regularization and the meson momentum cut-off A& in covariant regularization. As 
pointed out in Introduction, there are four parameters, the current quark mass mo, coupling 
constant G, quark momentum cut-off A/ and meson momentum cut-off A& to be fixed. By 
comparing with two observables m, = 139MeV, /„. = 92.4MeV and one reasonable empirical 
range of — 300MeV < 1/2 < qq > x / 3 < — 200MeV, we can not give fixed values of these four 
parameters. Here we regard the ratio z = A^/A^ as one free parameter. For each z, we 
can get a series of solutions from the above conditions. The meson cloud effect is now 
characterized by z, the larger z means more meson contributions. Specially, when z = 0, 
i.e., Ab = 0, it returns to the mean-field approximation. 

In the case of explicit chiral symmetry breaking mo 7^ 0, giving a z, for different m, we 
find a series of A^ and m from the above three equations, and then G, < qq > and other 
quantities could be calculated. In the case of chiral limit m$ = 0, giving a z, for different m, 
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we find Af and G from Goldberger-Treiman equation and the gap equation, then < qq > 
and other quantities are calculated. 

In both cases, the mesonic contributions to the quark self-energy m and to the pion 
polarization U n (k) are negative, and all the ratios of the mesonic corrections to the one- 
quark- loop contributions to m, H n (k) and physical quantities /„. and quark condensate 
increases with increasing z, and can reach 30-40% for values of parameters in the range 
0.3GeV < m < 0.5GeV and 1 < z < 1.5. These results are qualitatively the same as those 
H and JTJ . 



in 



Our numerical results are shown in Figs. 3-5, which concentrate on investigating the 
effect of explicit chiral symmetry breaking. In all figures, the solid lines correspond to the 
case of explicit chiral symmetry breaking m 7^ 0, and the dashed lines correspond to the 
case of chiral limit m = 0, and a, b, c, d, e correspond to z = 0, 0.5, 1, 1.5, 2, respectively. 

In Fig. ^, we plot the quark scalar condensate (—1/2 < qq > 1 ^ 3 ) as a function of 
constituent quark mass m for different values of z. This figure is different from the Fig. 6 
and Fig. 7 in [|T^]. Our numerical results show that: 

1) . For each curve, the quark condensate has a plateau, in which the quark condensate 
changes slowly with m, and the width of the plateau becomes more and more narrow with 
increasing z. For each curve, we define the plateau in the range of (—1/2 < qq > x / 3 ) < 
(-1/2 < qq >\2 n +0.0015) GeV. 

2) . For each z, the plateau of quark condensate is in the experimental range 0.2 GeV ~ 0.3 
GeV, and the plateau becomes higher and higher with increasing z, which shows that the 
meson contributions affect the quark condensate explicitly. 

3) . The solid and dashed lines of a are almost coincide, which means that the effect of 
current quark mass is invisible in the mean-field approximation. While beyong mean-field 
approximation, i.e., considering the meson contributions, the two curves in the case of b or 
c, d, e are separate in the region of the plateau. 
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FIG. 3. Quark condensate —1/2 < qq > 1//3 as a function of the constituent quark mass m for 



different values of z = A^/Aj. The solid lines and dashed lines correspond to mo / and ttiq = 



respectively. 



In Table 1, corresponding to the plateau region (—1/2 < qq > x / 3 ) < (—1/2 < qq > 
+0.0015)GeV, 





-| < QQ > 1/3 
(GeV) 


m(GeV) 


m (MeV) 


A/ (GeV) 


GA} 


z=0 


m = 


0.2096+0.0015 


0.47+0.1 




0.615+0.002 


4.8+0.9 


m 7^ 


8.9+0.30 


z=l 


m = 


0.2397+0.0015 


0.39+0.05 




0.725+0.01 


4.8+0.5 


m ^ 


0.2360+0.0015 


0.40+0.05 


7.78+0.20 


0.71+0.01 


4.7+0.5 


z=1.5 


m = 


0.2599+0.0015 


0.38+0.04 




0.801+0.008 


5.3+0.5 


mo^O 


0.2566+0.0015 


0.39+0.04 


7.29+0.20 


0.785+0.008 


5.2+0.5 



Table. 1. Quantities in the region of defined plateaus. 
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we list the range of consituent quark mass m, current quark mass mo, quark momentum 
cut-off Af and a dimensionless quantity GA^.This table shows that, in the mean-field ap- 
proximation, the values of constituent quark mass m and the current mass mo in the plateau 
is a little higher than the empirical values m ~ 1/3 proton mass and mo — 5 ~ 7MeV, and 
the quark condensate in the plateau is much lower than 0.25 GeV; however, these quantities 
in the plateaus at z = 1, 1.5 are more reasonable comparing with the empirical values. 

In Fig. we show the quark momentum cut-off Af as a function of the constituent 
quark mass m, 




FIG. 4. Quark momentum cut-off Af as a function of the constituent quark mass m for different 
values of z = A^/Af. The solid lines and dashed lines correspond to mo 7^ and mo = 
respectively. 

and in Fig. |5] we show the current quark mass mo as a function of the constituent quark 
mass m for different values of z. There is also a plateau for each curve. Although the 
ranges of constituent quark mass m corresponding to the Af plateau, mo plateau and quark 
condensate plateau are not exactly the same, they are mostly around m = 0.4 ± O.lGeV. 
This range m = 0.4±0.1GeV is determined dominantly by the gap equation in the mean-field 
approximation, the effects of meson cloud contributions are also shown clearly. Comparing 
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the plateaus at z = and z = 1.5, we see that Aj and (—1/2 < qq > 1 / 3 ) are all modified 
by the order of about 30%. But the corresponding plateaus at z = 2 modified too much. 
In fact, all these figures show that the plateaus of z = 2 are abnormal comparing that of 
z = 0.5, 1, 1.5. Based on the above analysis, it seems that z = 2 is not a reasonable choice, 
and z = 1 ~ 1.5 is more reasonable. 




0.001 I ' ' ' 1 ' ' ' 1 ' ' ' 1 ' ' ' 1 ' ' ' 1 ' ' ' 1 ' ' ' — I 
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FIG. 5. Current quark mass mo as a function of the constituent quark mass m for different values 
of z = Afj/Af. The solid lines and dashed lines correspond to uiq ^ and tuq = respectively. 

In the region m = 0.4 ±0.1 GeV, all these three figures Figs. |3|-|5| show that: 

1) In the mean-field approximation, the solid line and dashed line of z = coincide, which 
means that the effects of the current quark mass are invisible in this case. 

2) When beyond mean-field approximation z > 0, i.e, considering meson contributions, the 
effects of current quark mass are explicitly reflected by the separation of the solid and dashed 
lines. 

As we pointed out in Introduction, in NJL model all information is calculated from con- 
stituent quark propagator, and the current quark mass mo induces the massive pion m n . In 
the mean-field approximation, there is no feedback of meson modes to quark propagator, so 
the effects of m or m n can not be reflected. Only beyond mean-field approximation, when 
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the back interaction of meson mode to quark propagator are considered, can the effects of 
m or m n be reflected on the quantities calculated from the fundamental element, quark 
propagator. 

V. CONCLUSIONS AND SUMMARY 

We establish a general scheme of SU(2) NJL model including the current quark mass 
explicitly and considering the meson cloud contributions. To ensure the quark self-energy 
expanded in the proper order of 1/N C , we discuss in details that the internal meson propaga- 
tor must be in O(l/N ) order. And we also develop a unified way to calculate all Feynman 
diagrams in the NJL model by only expanding the constituent quark propagator in external 
momentum k. Our numerical results show that: a) when z — 1 ~ 1.5, the NJL parameters 
are more reasonable than that in the mean-field approximation; b) the effect of current 
quark mass can be reflected only beyond mean-field approximation, and it is characterized 
explicitly by the separation of the solid and dashed lines in the range of the parameters 
corresponding to m between 0.3GeV ~ 0.5GeV. 
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APPENDIX: Definition of Integral Functions 
and Corresponding Pauli-Villars Expressions 



The fundamental functions are defined as following: 

d 4 p 1 



K(q) 



L{q) 



M{q) 



(27r)V-m 2 ' 
d 4 p 1 



(2vr) 4 (p 2 


— m 2 )((p + q) 2 — m 2 ) ' 


• d 4 p 


1 


(2tt) 4 {p 2 


— m 2 ) 2 ((p + q) 2 — m 2 ) 


d 4 p 


1 


(2tt) 4 {p 2 - 


- m 2 ) 2 ((p + q) 2 — m 2 ) 2 


<• d 4 p 


1 



(27r) 4 (p 2 — m 2 ) 3 ((p + q) 2 — m 2 ) ' 
And their corresponding Pauli-Villars expressions are: 

M 2 



1 PV 



' £<7 a M a 2 lli% 

(An) 2 ^ m 2 



I(q)pv = j^y 2 E C«(- In ^ + 2(1 - ^1 + y" 1 ) ln(^ + V /T + ^), 



1 ln( v ^+ v ^Ty:) 



Wry ~E^(- 8Ma4 J 1 + 2/a) )(l - -^^ln(VK + ^)), 



M\/^ + a/ 1 + Va) 



(4tt)^ aV 16M 4 y a (l + y a ) yV ^ ^y fl (l + y a ) 

where C a = [1, 1, -2], a a = [0, 2, 1], M 2 = m 2 (l + a a x) for a=0,l,2, and y a = -g 2 /( 4 - 
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